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In a superconductor coupled with a ferromagnetic metal, spin and charge imbalances can 
be induced by injecting spin-polarized electron current from the ferromagnetic metal. We 
theoretically study a nonequilibrium distribution of quasiparticles in the presence of spin 
and charge imbalances. We show that four distribution functions are needed to characterize 
such a nonequilibrium situation, and derive a set of linearized Boltzmann equations for 
them by extending the argument by Schmid and Schon based on the quasiclassical Green's 
function method. Using the Boltzmann equations, we analyze the spin imbalance in a thin 
superconducting wire weakly coupled with a ferromagnetic electrode. The spin imbalance 
induces a shift 5/i (—Sfi) of the chemical potential for up-spin (down-spin) quasiparticles. 
We discuss how <5/i is relaxed by spin-orbit impurity scattering. 

KEYWORDS: nonequilibrium superconductivity, spin imbalance, charge imbalance, spin-orbit 
interaction 

1. Introduction 

Electron current injected into a superconductor produces a nonequilibrium distribution of 
quasiparticles. 1-9 A number difference between electrons and holes arises in such a nonequi- 
librium situation. 4,5 This is called charge imbalance. The charge imbalance induces an ex- 
cess quasiparticle current, which results in a potential difference between pairs and quasi- 
particles. 1 ' 4 ' 5 Although early experiments on the charge imbalance have focused on spin- 
independent phenomena, we can now study spin-related ones by injecting spin-polarized elec- 
tron current into superconductors by using a ferromagnetic metal as an injection electrode. 10 
If injected electron current is spin polarized, there arises spin imbalance between up- and 
down-spin quasiparticles in addition to the charge imbalance. In this case, a potential dif- 
ference between up- and down-spin quasiparticles arises in addition to a spin-independent 
potential difference due to the charge imbalance. The spin imbalance in mesoscopic supercon- 
ductors has recently attracted much attention both experimentally 11-13 and theoretically. 14,15 
The main attention is focused on its relaxation due to spin-flip processes. 

The most general framework to understand the charge imbalance in superconductors is 
presented by Schmid and Schon. 7 Based on kinetic equations for the quasiclassical Green's 
functions, they showed that nonequilibrium quasiparticle distributions are described by two 
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distribution functions f^(r,e) and fi,(r,e), where e represents quasiparticle energy measured 
from the chemical potential [i at equilibrium, and derived a set of linearized Boltzmann equa- 
tions for them. It has turned out that /t and /l describe the charge imbalance and the related 
excess quasiparticle energy, respectively, and that they are coupled with respectively trans- 
verse and longitudinal variations of the pair potential. By using the Boltzmann equations, 
we can describe the charge imbalance including effects of inelastic phonon scattering. How- 
ever, their framework is restricted to spin-independent phenomena, and cannot apply to the 
case where the spin imbalance plays a role. To fully understand nonequilibrium quasiparticle 
distributions in the presence of a spin-polarized current injection, we need a framework by 
which we can describe both the spin and charge imbalances in a unified manner including 
effects of phonon scattering. Although a few theoretical treatments for the spin imbalance in 
superconductors have been reported so far, 14,15 they do not satisfy all our requirements. 

In this paper, we study a nonequilibrium quasiparticle distribution in the presence of 
the spin and charge imbalances. To describe such a nonequilibrium situation, we adopt the 
kinetic equation approach by Schmid and Schon. We introduce four distribution functions 
/l-j /t+ an d /t- to characterize quasiparticle distributions. It is shown that /l+ and /t- 
represent the spin and charge imbalances, respectively, and /l- and /t+ represent the excess 
quasiparticle energy and the energy imbalance between up- and down-spin quasiparticles, 
respectively. The suppression of the pair potential due to a current injection is determined by 
/l- . We derive a set of linearized Boltzmann equations for the distribution functions in steady 
states assuming that the spin imbalance is relaxed by spin-orbit impurity scattering. We show 
that the four distribution functions are decoupled with each other in the resulting Boltzmann 
equations. This indicates that we can separately consider the spin and charge imbalances in 
steady states. As an application of the Boltzmann equations, we treat the spin imbalance in 
a quasi-one-dimensional superconducting wire weakly coupled with a ferromagnetic electrode 
through a tunnel junction. The spin imbalance induces a shift 5fi {—SfJ,) of the chemical 
potential for up-spin (down-spin) quasiparticles. We discuss spatial decay of 5/j, due to spin- 
orbit impurity scattering. At high temperatures, the decay of Sfjb(x) obeys the exponential 
law e - l x-x °l/ As , where A s is the spin-diffusion length and \x — xq\ represents the distance 
from the injection point xq. This is in agreement with the previous result. 15 However, at low 
temperatures, we observe deviations from the exponential law near the injection point. 

In the next section, we consider a superconductor coupled with a ferromagnetic metal 
and introduce the quasiclassical Green's functions for the superconductor in the Keldysh 
formalism. The kinetic equations for them are derived in the presence of spin-orbit impurity 
scattering. In §3, we introduce four distribution functions to describe both the spin and charge 
imbalances, and derive a set of Boltzmann equations. We clarify the meaning of each distri- 
bution function. In §4, we analyze the spin-imbalance relaxation in a thin superconducting 
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wire based on the resulting Boltzmann equations. Section 5 is devoted to a short summary. 
We set h = &b = 1 throughout this paper. 

2. Kinetic Equations for Green's Functions 

We consider a superconductor coupled with a ferromagnetic metal through a point-like 
tunnel junction. As a model for the superconductor, we adopt an electron gas interacting with 
phonons. Phonons are assumed to be described by the Debye model with the sound velocity 
Co- We assume that electrons experience normal impurity scattering and spin-orbit impurity 
scattering. Let ij) a (x) with x = (r,t) be the electron field operator in the superconductor. We 
introduce = t (ip^(x),t/j^(x)) and define the following Green's functions 16 

G K (x,x') = -if z ([*(x),& (*')]_), (1) 

G R (x,x>) = -if z S(t- 0< [*(*), *V)]+>, (2) 

G A (x, x>) = +if z e(t' - * f (*')] +>, (3) 

where @(t) is Heavisid's step function. Here and hereafter, fj (i = x, y, z) represents the Pauli 
matrix. We use the Keldysh representation 

r( i\ ( G R (x,x>) G K (x,x>) \ 

G(x,x ) = , (4) 

V G A (x,x') J 

and define its Fourier transform as 

G(r,p,t,t') = |d 3 S e-^G(r + |,t,r-|,t'). (5) 
Integrating this over ^ = p 2 /(2m) — /j,, we obtain the quasiclassical Green's function 17 

G(r,p,t,t') = ^ j dHG(r,p,t,t'), (6) 
where p = p/\p\- It has been shown that the quasiclassical Green's function obeys 18 ' 19 
v F p ■ VG(r, p, t, t') + r z d t G{r , p, t, t') + fyG(r, p, t, 

/oo 
dti(E(r,p,t,*i)G(r,p,ti,0-G(r,p,t,*i)E(r,p,*i,0) =0, (7) 
-oo 

where T i = diag(fj, fj) and S(r,p,t,i') represents the self-energy part. The self-energy part 
is decomposed into 

E(r , p, t, t') = S imp (r, p, t, t') + ^(r , p, t, t') + S ph (r , p, t, t') + S inj (r, p, t, t'), (8) 

where S imp , S so and S ph represent the contributions from normal impurity scattering, spin- 
orbit impurity scattering and electron-phonon interaction, respectively. It should be noted 
that S ph contains the pair potential. The last term S in j represents the spin-polarized current 
injection from the ferromagnetic metal. We simplify eq. (7) by following the argument by 
Usadel. 20 We employ an approximation 

G(r,p,t,t') = G(r,t,t')+p-G 1 (r,t,t'), (9) 
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where G_i(r,t,t') is much smaller then G(r,t,t'). After some manipulations, we obtain 
DV • (G(r, e, t)VG(r, e, t)) + ie [r 2 , G(r, e, t)]_ 

- ld t \r z ,G(r,e,t)} + - i[Ur,e,t),G(r,e,t)]_ = 0, (10) 

where 

G(r,e,*) = Jdse Us G(r,t + ^t-^). (11) 

In the following, we restrict our consideration to steady states, and neglect the ^-dependence 
of all the functions. 

We evaluate the elements of the self-energy part. We first consider S imp and £g . The 
Hamiltonian H imp for normal impurity scattering is written as 

ffimp = E / d3r VtoVW WVv(r). (12) 
We assume for simplicity that the impurity potential is given by 

Vi m p(r) = u imp £ 5(r - Rj), (13) 

where Rj indicates the position of the jth impurity. The Hamiltonian H so for spin-orbit 
impurity scattering is written as 

H so = Y.j dV ^l(r)U^'(r)^,(r). (14) 

<t,<t' 

The spin-dependent potential is 

U°f\r) = ■ (W so (r) x V) , (15) 

where V so (r) is given by eq. (13) with the replacement of Ui mp — ► u so /kp. In this case, S imp 
and S so are given by 

S imp (r, e ) = --^G(r,e), (16) 

^ T imp 

S so (r, e) = -77^G(r, e) + J-T^d{r, e)r y . (17) 

D^so "''"so 



The relaxation times are given by 

1 



= 27m imp 7Vs(0)|u imp r, (18) 

''"imp 

— = 2vm imp iVs(0)K o | 2 (Pxp)i. ( 19 ) 

''"so 

i=x,y,z 

where n\ mp and iVg(O) represent the impurity concentration and the density of states in the 
normal state, respectively. Note that the second term of S so with G_' represents spin-flip 
processes. The Green's function G'(r,e) is the Fourier transform of G'(x,x'), where G_'(x,x') 
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is defined by eqs. (l)-(4) with the replacement ^>(x) — > ^' (x) = ipi(x)). We can show 

that its elements G x (x,x') (X = R,A,K) satisfy 

G' x (x,x') = -f z t [G x (x',x)]f z , (20) 

where R = A, A = R and K = K. These relation will be used later. We obtain 

S imp (r, e) + ^(r, e) = ~G(r, e) - {-r y d(r, e)r y - G(r, e)) (21) 

with t- 1 =Tr^ + T-K 

We turn to S in j. Let us consider the case where a ferromagnetic metal is coupled to the 
superconductor through a point-like tunnel junction at Vq. The corresponding Hamiltonian 
#inj is given by 

#inj = E / d ' VdV ( T ( r ' r'WUr)Mr')e- ieVt + h.c.) , (22) 

where 4> a and V represent the electron field operator in the ferromagnetic metal and an applied 
bias voltage, respectively, and we assume that T(r,r r ) = T5{r — ro)6(r' — r'o)- With this 
model, we obtain 

£? Q} (r,e) = -i\T\ 2 5(r-r )a x i (e) (23) 

with 

where A?p-f(0) (iVpjXO)) represents the density of states for up-spin (down-spin) electrons in 
the ferromagnetic metal. 

The phonon self-energy part has been presented in ref 18. Note that T,^ h and Sp h contain 
the pair potential, so we separate it out as 

t^ A (r, e) - ^(r, e) - if x (r)A(r) (26) 

with 

f x (r) = cos x(r)f y + sinx(r)f x , (27) 

where x( r ) an d A(r) represent the phase and amplitude of the pair potential, respectively. 
Assuming that phonons are in thermal equilibrium, we obtain 

s2rV,6) = -i J deV ph (6, e '){coth(^)G^(r, e ')T^(r,e')}, (28) 



S£(r,e) = -i 



| d e V ph (e,e'){coth(^)G^(r,eO- (G R (r,e')-G A (r,e'))}, (29) 
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where 

^ P h(6, 6') = • j^yV ~ e ) 2 sign(e' - e). (30) 

Here, Co and g represent the sound velocity and the coupling constant for the electron-phonon 
interaction, respectively. 

Using eq. (10) and the elements of the self-energy part, we obtain the kinetic equations 
for the Green's functions, 

DV ■ (6 R ' A VG R ' A } + ie[f 2 , G R,A ] _-A[f x , G R > A ] _ + f R > A - I R > A - P R > A = 0, (31) 

DV ■ (G R VG K + G K VG A ) + ie [t z , G K ] _ - A [f x , G K ] _ + f K - I K - P K = 0, (32) 
where 

t R,A = 1 ^ y G'R^r y G R ' A - G R ' A T y 6' R ' A Ty) , (33) 



f K = -L (f y G' R r y G K + f y G' K f y G A - G R f y G' K f y - G K r y G' A f y ) , (34) 
3t so 

i R ' A = i(t R ^ A G R ' A - G R ' A ± R ^ A ) , (35) 

I K = ^%G K + t^G A - G R t% - G K t A ) , (36) 

and P x is obtained by the replacement of — > in I x . The arguments r and e are 
suppressed in the above equations. Here, T x represent the influence of spin-flip processes due 
to spin-orbit impurity scattering, and I x and P x describe inelastic phonon scattering and a 
spin-polarized current injection, respectively. If we neglect T x and set A^f| (0) = A^p|(0) in P x , 
our argument is reduced to the previous one presented by Schmid and Schon. 7 Quasiparticle 
distribution functions are contained in G K and G K . 

3. Boltzmann Equations 

Based on the kinetic equations presented in the previous section, we start to derive a set 
of linearized Boltzmann equations which describes a nonequilibrium quasiparticle distribution 
in superconductors. In terms of the spectral functions Aq, A^2, R\ and R2, we approximately 
express G R and G A as 7 

G R ' A {r,e) = (±N 1 (r,e)+iR 1 (r,e))f z + (N 2 {r,e)±iR 2 {r,e))r x {r). (37) 

The spectral functions satisfy 

N h2 (r,e) = N h2 (r,-e), (38) 

Ri, 2 (r,e) = -R h2 {r,-e). (39) 
We rewrite f y G x (r, e)f y in terms of G x (r, e). From eq. (20), we obtain 

G' x (r,e) = -f z t [G x (r,-e)]f z . (40) 
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Combining eqs. (37) and (40) with eqs. (38) and (39), we obtain T y G R (r, e)f y = —G R (r, e) and 
T y G ' A (r , e)f y = —G A (r,e). Note that T R > A = if these results are substituted into eq. (33). 
This indicates that spin-flip processes due to spin-orbit impurity scattering do not influence 
on the spectral properties. 21 Only quasiparticle distribution functions contained in G K are 
affected by spin-flip processes. We adopt the following expression 19 

G K (r, e) = G R (r, e)h(r, e) - h(r, e)G A (r, e) (41) 

with 

h(r,e) = h 1 (r,e) + h 2 (r,e)t z . (42) 

A nonequilibrium quasiparticle distribution is described by hi and h 2 . From eq. (40), we 
obtain 

T y G' K (r, e)t y = G R (r, e)f y h{r, -e)t y - t y h(r, -e)f y G A (r, e). (43) 
By using eqs. (41) and (43), we simplify the expression of as 

r K (r, e) = ^T (G R (r, e)5h(r, e)G A (r, e) - 5h(r, e)) (44) 

with 

6h(r, e) = (h^r, e) + h^r, -e)) + (h 2 (r, e) - h 2 (r, -e)) t z . (45) 
At equilibrium, hi^{r,e) is reduced to hi = tanh(e/(2T)) and h 2 = 0. It is convenient to set 



/»i(r,€)=tanh( k — J -2/ L (r,e), (46) 

h 2 (r,e) = -2f T (r,e). (47) 

The quasiparticle distribution functions are written as 7 

Mr, e) = / FD (e, n) + h(r, e) + / T (r, e), (48) 

ft(r, e) = /fd(c, /i) - h(r, -e) + / T (r, -e), (49) 

where /fd(£,A*) is the Fermi-Dirac distribution function. It should be emphasized that 
fh(f,e) = —fi,(r,—e) and /x(r, e) = /^(r,— e) are implicitly assumed in ref. 7. These re- 
lations straightforwardly result in /|(r,e) = /j.(r, e). Thus, the framework by Schmid and 
Schon is restricted to spin-independent phenomena. We do not accept the symmetry relations 
to enable us to consider the spin imbalance. 

We obtain Boltzmann equations based on eqs. (31) and (44). Substituting eqs. (41) and 
(43) into eq. (32), we derive equations for e) and /t(t", e) with the help of eq. (31). 
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Details of the derivation are described in ref. 16. We obtain 

D • V(iV 2 (r, e) - R 2 2 (r, e)) V/ L (r, e) - ^- (iV 2 (r, e) - R 2 2 (r, e)) (/ L (r, e) + / L (r, -e)) 

+ / L (r,e,{/ L }) +P L (r,e) = 0, (50) 
L> • V(AT 1 2 (r, e) + iV|(r, e)) V/ T (r, e) - ^-{N 2 (r, e) + iVf (r, e)) (/ T (r, e) - fr(r, -e)) 

Vr /T ( r ' £ ) + M r > e > < /t K> + PT ( r ' e ) = °' ( 51 ) 

T<xmv\fi) 

where r CO nv(e) represents the conversion time for charge imbalance. 22 The collision integrals 
7l and It due to inelastic phonon scattering are expressed as 

Il,t{v, €,{/}) =-2 y deV ph (e,e / )M LiT (r,e,e / ) 

cosh 2 (#) /(r, e) - cosh 2 (^)/(r, e') 



sinh(^y^) cosh (27) cosh(^) 



(52) 



where 



M L (r, e, e') = ^(r, e)Wi(r, e') - R 2 (r, e)R 2 (r, e'), (53) 
M T (r,e,e') = N^r^^r^') + N 2 (r,e)N 2 (r,e'). (54) 
The injection terms Pl and Pt are given by 

,e-eV, 



P L (r,e) = ||T| 2 5(r-r )iV 1 (r, e ){iV FT (0)(tanh(^)-tanh(- ^ 

+ iV Fi( 0)(tanh(^) -tanh(i^))}, (55) 
P T (r, e) = ||T| 2 <5(r - r )7Vi(r, e){iV FT (0) (tanh(^) - tanh( 



e-eV 
2T 



-iV Fi (0)(tanh(^) -tanh(i±^))}. (56) 



To make the equations much simpler, we introduce the four distribution functions, 

/ L+ (r, e ) = i(/ L (r,e) + / L (r,-e)), (57) 

/L-(r,e) = I(/ L (r,e)-/ L (r,- e )), (58) 

/ T+ (r, e) = \ (Mr, e) + / T (r, -e)) , (59) 

/T-(r,e) = i(/ T (r,e)-/ T (r,-6)). (60) 
We observe that they satisfy 

/L,T+(r, -e) = / L ,T+(r, e), (61) 

/ Li T-(r,-e) = -/ L , T -(r,e). (62) 
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Noting eqs. (30), (38) and (39), we obtain a set of Boltzmann equations for /l± and /t± as 
D-V(N 2 (r,e)-R 2 (r,e))Vf L+ (r,e) 

- — (JV?(r, e) - R 2 2 (r, e))/ L+ (r, e) + J L (r, e, {/ L+ }) + P L+ (r, e) = 0, (63) 
D ■ V(N 2 (r, e) - R 2 2 (r, e)) V/ L -(r, e) + J L (r, e, {/ L _}) + P L _(r, e) = 0, (64) 

D • V(AT 1 2 ( r , e) + N 2 (r, e)) V/ T+ (r, e) ^/r+(r, e) + J T (r, e, {/ T +}) + Pr+(r, e) = 0, 

(65) 

D • V(AT 1 2 ( r , e) + iV 2 (r, e)) V/ T -(r, e) ^/r-(r, e) 

7"conv V^J 

- — (JV?(r,e) + iV|(r,e))/ T _(r,e) + / T (r, e, {/ T _}) +Pr_(r,e) = 0, (66) 

Tsf 

where 1 = (4/3)r S o 1 and 

PL±(r,e) = ^(P L (r,e)±P L (r,-e)), (67) 

**r±(r, e) = \ (Pr(r, e) ± Pr(r, -e)) . (68) 
Using the expression of the tunnel resistance, 

^ 1 = 2^e 2 iVs(0)(iV FT (0) + N n (0))\T\ 2 , (69) 
we can rewrite the injection terms as 



P -(-» = S« a (- - "«>M^) - tanh(l^)), (72) 

Pt -<"' £ » = - -»)( 2 ^^) - tanh(i±^) - tanh(^)), (73) 

where the spin polarization P s is defined by 

= iV FT (0)-iV Fi (0) 
s iV FT (0) +N Fl (0)' 1 j 

Equations (63)- (66) and eqs. (70)- (73) are the central result of this paper. We here clarify 
the meaning of each distribution function. The distribution functions f a (r,e) are expressed 
as 

Mr, e) = / FD (e, /*) + f L+ (r, e) + / L _(r, e) + / T+ (r, e) + / T -(r, e), (75) 
Mr,e) = / FD (e,/x) - / L+ (r,e) + / L -(r,e) + / T+ (r,e) - / T -(r,e). (76) 
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Let S(r) and Q(r) be the spin and charge imbalances, respectively. Noting that Ni(r,e) is 
the normalized local density of states in the superconductor, we can express 

/oo 
deAq(r,e)(/ T (r,e)-/;(r,e)), (77) 
-oo 

/oo 
de Aq(r, e) (/ T (r, e) + f^r, e) - 2/ FD (e, //)) . (78) 
-oo 

By using eqs. (75) and (76), we obtain 

/•oo 

S(r) = 4iV s (0) / de N^r, e)/ L+ (r, e), (79) 
■/ o 

/•OO 

Q(r) = 4AT S (0) / deiVi(r,e)/r+(r,e). (80) 
■/ o 

Thus, /l+ and /t+ describe the spin and charge imbalances, respectively. Other two distribu- 
tion functions are related to quasiparticle energies. Let Eq and Eg be the excess quasiparticle 
energy and the energy imbalance between up- and down-spin quasiparticles, respectively. They 
are given by 

/oo 
de N^r, e)e (/ T (r, e) + f^r, e) - 2/ FD (e, fi)) , (81) 
-oo 

/oo 
deJVi(r,e)e(/ T (r,e)-/i(r,e)). (82) 
-oo 

By using eqs. (75) and (76), we obtain 

/•oo 

E Q (r) = 4iV s (0) / de Aq(r, e)e/ L _(r, e), (83) 
J o 

poo 

E s (r) = 4N s (0) deN 1 (r,e)ef T -(r,e). (84) 
Jo 

Thus, /l- and /x_ describe the excess quasiparticle energy and the energy imbalance be- 
tween up- and down-spin quasiparticles, respectively. The two distribution functions /l+ and 
/t- characterize the spin imbalance, while other two distribution functions /l- and /t+ 
characterize the charge imbalance. The former two have not been discussed in literatures. 

We approximately obtain the spectral functions. The presence of spin-orbit scattering 
does not result in any changes of the spectral properties of superconductors. According to the 
approximation adopted in eq. (37), we obtain 

A(r) = N s (0)g 2 J™deR 2 (r,e) (tanh(^) - 2/ L _(r, e)) . (85) 

This indicates that the suppression of the pair potential due to a current injection is governed 
by /l-. In contrast, the variation of the phase x 1S related to /x+ although this point is out 
of our scope. We assume that A and x spatially vary much slower than /l± and /x±- Thus, 
we approximate that X7G R ' A = 0. Furthermore, we neglect the phonon self-energy in deriving 
G R ' A . The local density of states vanishes for |e| < A(r) in this case, so we consider the energy 
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region of |e| > A(r) in the following. After these simplifications, we obtain 



Aq(r,e) 



e 



(86) 



R2{r,e) 



^2 _ A 2( r ) 

sign(e)A(r) 



(87) 



^ e 2 _ A 2( r ) 



and N 2 (r,e) = for |e| > A(r). 

In the presence of the spin and/or charge imbalances, the distribution functions f a (r,e) 
deviate from the equilibrium ones. To characterize their deviations, we introduce the spin- 
dependent chemical potential /v(r) = M + <W(f) for quasiparticles. To define <5/v(r), we 
assume that a fictitious electrode is weakly coupled to a superconductor at r through a point- 
like tunnel junction. In terms of a bias voltage Vfi c , the spin-dependent tunneling current 
Ia{r, Vfi c ) is given by 



The tunneling current I a (r, Vfi c ) vanishes if eVg c is equal to Sfx a (r). Thus, 6[i a (r) satisfies the 
following equation 



We use the above equation as the definition of S/j, a (r). We observe that <5/x-f(r) = <5/ij_(r) if 
the spin imbalance is absent (i.e., /L+(r,e) = 0), while 8fx^{r) = — <5/^(r) in the absence of 
the charge imbalance (i.e., /x+(r, e) = 0). 

4. Spin-Imbalance Relaxation 

In this section, we study the behavior of spin-imbalance relaxation based on the Boltzmann 
equation for Let us consider a thin superconducting wire coupled with a ferromagnetic 
electrode through a tunnel junction. We assume that /l_ is very small everywhere in the 
superconductor, and set A(r) = A. We neglect the charge imbalance for simplicity, and focus 
on a shift of the spin-dependent chemical potential. If the cross-sectional area A of the wire 
is small enough, we are allowed to consider a one-dimensional problem of /l+(x, e). As noted 
just below eq. (89), we observe that Sfi(x) = 8^{x) = —5^{x) in this case. The Boltzmann 
equation for /l+ is reduced to 





(88) 





(89) 



Ddlf-L+(x,e) - —f L+ (x,e) + J L (x, e, {/l+}) +-FL+(x,e) = 



(90) 
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where 

I L (x,e,{f}) = -2 j deV ph (e,e') (N^N^e') - R 2 (e)R 2 (e')) 
cosh 2 (27) f(x, e) - cosh 2 (^) f(x, e') 



sinh(^^) cosh^) cosh(^) 



(91) 



We define the energy relaxation time T cnc (e) as 

1 (' coshf— ) 

—— = 2 / deV ph (e,e') (N 1 (e)N 1 (e') - R 2 (e)R 2 (e')) . ^ > ( e , y (93) 
^ene(e) J smh(^) cosh(^,J 

in terms of which the first term in the collision integral /l is rewritten as — /L+(r, e)/r cnc (e). 

Note that the energy of an injected quasiparticle is within A < |e| < eV + T. The behavior of 

the spin- imbalance relaxation depends on whether r s f is longer or shorter than r ene (e) in this 

energy range. 

We first consider the high-temperature regime in which eV <^T and T ene (e) is much shorter 
than r s f for A < |e| < T. In this case, the energy dependence of /l+ is mainly determined by 
the collision-integral term. Except near the injection point (i.e., x = xo), we can approximate 
fL + (x,e) oc 1/ cosh 2 (e/(2T)). Since eq. (89) is simplified to 

we observe that 

f L+ (x,e) = (95) 

J +K ' 4Tcosh 2 (^) V ; 

We determine 5fi(x) based on eq. (90). If we approximately set iVi(e) = 1 in P\ J+ (x,e), we 
obtain 

5fi(x) = 5jiQ e A = ; (96) 



where A s = \JDt s { is the spin-diffusion length and 

2X s P s eV 

Sf *° = 8e^N s (0)R t AD- (97) 
The approximation Aq(e) — > 1 results in an under-estimation of 5[1q. It should be noted that 
the distribution function /|,j.(x, e) = /fd(c, jj) ± fh+{e, n) is expressed by the Fermi-Dirac 
distribution function with a shifted chemical potential as e) pa /fd(c, fi±5/j,(x)) in this 

case. 

Next, we consider the low-temperature regime in which T ene (e) is much longer than r s f for 
A < |e| < eF. Thus, we can neglect the collision-integral term in eq. (90). Solving eq. (90), 
we obtain 

f L+ (x,e) = f (e)e- t ^ ni (98) 
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Fig. 1. S/i/A for T/A = 0.04 and 0.08 as a function of the normalized distance y = \x — x \/\ s from 
the injection point. The solid and dotted lines correspond to eV/A — 1.2 and 1.6, respectively. 

with 

/o ( e ) = 9A7 A £ S „ (e) ( tanh(iif^) - tanhf^^ . (99) 

,/uv 7 16e 2 iV s (0)i?iAD w V 2T 2T / 

In this case, we cannot express the distribution function f a {x, e) by the Fermi-Dirac distribu- 
tion with a shifted chemical potential. Substituting eq. (98) into eq. (89), we obtain 

jTd £ Wl(£ )(tanh(i±|jM) - «anh(^|M) - 4AM e"^) = 0. (100) 

We numerically solve eq. (100) and obtain 5fi(x)/A as a function of the normalized distance 
y = \x — xo|/As from the injection point xq at eV/A = 1.2 and 1.6 for T/A = 0.04, 0.08 
and 0.16. The following parameters are adopted: Rt = 2 kS7, iVg(0) = 1.2 x 10 22 eV _1 cm~ 3 , 
D = 5.3 x 10 9 ^m 2 /s, r s f = 200 ps and A = 50 x 250 nm 2 . The values of D and r s f result in 
A s = 1.03 fim. Figure 1 shows that the bias- voltage dependence of 5fi is very weak, in contrast 
to the high-temperature regime where 5/x oc eV . This clearly indicates a nonlinear nature of 
eq. (100). From Fig. 2, we observe that the decay of 5fi(x) can be fitted by the exponential 
law e - l x-x °l/* s in the case of T/A = 0.16. However, except for this case, a deviation from 
the exponential law appears near the injection point (i.e., y = 0). We also observe that the 
asymptotic behavior of 5^,{x) far from the injection point is still governed by the exponential 
law. It should be noted here that the anomalous slow decay of 5fi observed near the injection 
point is partly attributed to the divergence of Aq(e) at the gap edge. Since the divergence is 
smeared by a gap anisotropy, the anomalous behavior may be weakened in actual cases. 

Yamashita et al. 15 have studied the spin-imbalance relaxation in superconductors by as- 
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Fig. 2. Sfi/A at eV/A = 1.2 as a function of the normalized distance y = \x — x |/A s from the 
injection point. The solid, dotted and dashed lines correspond to T/A = 0.04, 0.08 and 0.16, 
respectively. 



suming that the quasiparticle distribution function can be expressed in the form of the Fermi- 
Dirac distribution function with a shifted chemical potential n a (x) = /U±<5/x(x), and concluded 
that 5[i a (x) decays exponentially on the length scale of A s . Our argument indicates that the 
assumption employed in ref. 15 can be justified only in the high-temperature regime. Indeed, 
near the injection point at low temperatures, we have found deviations from the exponential 
law. 

We have neglected the charge imbalance in the above argument. The nonlinear equation, 
eq. (89), is reduced to 

^ u (x) = 4Tcosh 2 (^)(±/ L+ (x,e) + / T+ (x,6)) (101) 

in the high-temperature regime. Thus, an additional spin-independent correction is simply 
added to the chemical potential if we take the charge imbalance into account. However, such 
a simple treatment cannot be applied to the low-temperature regime because we must solve 
eq. (89) in its present form to obtain 5ji a . Thus, the spin- imbalance and charge-imbalance 
corrections to 5[x a are not necessarily additive. 

5. Summary 

We have studied a nonequilibrium distribution of quasiparticles in the presence of both 
spin and charge imbalances. By extending the kinetic equation approach by Schmid and Schon 
based on the quasiclassical Green's function method, we have presented a set of linearized 
Boltzmann equations for distribution functions fh±(r,e) and fi±{r,e) in steady states. It is 
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shown that /l+ and /t- represent the spin and charge imbalances, respectively, and /l- and 
/t+ represent the excess quasiparticle energy and the energy imbalance between up- and down- 
spin quasiparticles, respectively. It is also shown that the suppression of the pair potential 
due to a current injection is governed by /l-- These distribution functions are decoupled with 
each other in the Boltzmann equations. This allows us to separately consider the spin and 
charge imbalances in steady states. 

As an application of the Boltzmann equations, we have considered the relaxation of spin 
imbalance in a quasi-one-dimensional superconducting wire weakly coupled with a ferromag- 
netic electrode. The spin imbalance induces a shift 5fi {—5 fx) of the chemical potential for 
up-spin (down-spin) quasiparticles. We have analyzed spatial decay of Sfi due to spin-orbit 
impurity scattering. We have shown that at high temperatures, 5fx decays exponentially on 
the length scale of A s , where A s is the spin diffusion length. However, at low temperatures, we 
have observed deviations from the exponential law near the injection point. 
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